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ABSTRACT: The implication of primordial magnetic-field-induced structure formation for the HI 
signal from the epoch of reionization is studied. Using semi-analytic models, we compute both the 
density and ionization inhomogeneities in this scenario. We show that: (a) The global HI signal 
can only be seen in emission, unlike in the standard ACDM models, (b) the density perturbations 
induced by primordial fields, leave distinctive signatures of the magnetic field Jeans' length on the 
HI two-point correlation function, (c) the length scale of ionization inhomogeneities is ^ 1 Mpc. 
We find that the peak expected signal (two-point correlation function) is ~ 10~ 4 K 2 in the range 
of scales 0.5-3 Mpc for magnetic field strength in the range 5 x 10~ 10 -3 x 10~ 9 G. We also 
discuss the detectability of the HI signal. The angular resolution of the on-going and planned radio 
interferometers allows one to probe only the largest magnetic field strengths that we consider. They 
have the sensitivity to detect the magnetic field-induced features. We show that the future SKA has 
both the angular resolution and the sensitivity to detect the magnetic field-induced signal in the 
entire range of magnetic field values we consider, in an integration time of one week. 
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1. Introduction 

Understanding the epoch of reionization is one of the outstanding challenges of modern cosmol- 
ogy. In recent years, a partial understanding of this important stage in the history of the universe 
has emerged. Gunn-Peterson tests have revealed the existence of neutral hydrogen (HI) at red- 
shifts z ^ 5.7 ([Q, E$L ^ ||]). The detection of temperature-polarization cross-correlation and the 
polarization-polarization correlation at large angular scales {i ^ 10) in the WMAP data have given 
firm evidence that the universe reionized around z ~ 10 (e.g. 

The most direct approach to studying this epoch is to attempt to detect the redshifted 21 -cm 
line of the neutral hydrogen during the period when the universe makes a transition from being fully 
neutral to fully ionized. There is substantial on-going effort in that direction and this also remains 
one of the primary goals of upcoming and future radio interferometers (e.g. MWA, for a detail de- 
scription see LOFAR , and SKA 2 ). Primordial tangled magnetic fields can substantially alter 
the scenario of structure formation of the universe hereafter Paper I and references therein). It 
has been shown that magnetic field-induced structure formation can cause early reionization (Pa- 
per I, [pt]) and might leave detectable signatures in the HI signal from the epoch of reionization 



cqiotEjji). 

In this paper, we study the implications of the existence of the primordial magnetic fields for 
possible studies of the reionization epoch. We compute both the global HI signal and its fluctuating 
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component. In particular, we study in detail the nature of ionization inhomogeneities in the pres- 
ence of magnetic fields, as their precise understanding is needed to correctly set the scale of the 
effect. The detailed quantification of this aspect constitutes the main difference between our present 
work and earlier studies. We also discuss in detail the detectability of the magnetic field-induced 
HI signal. 

In §2, we briefly review the salient features of the structure formation process in the presence 
of primordial tangled magnetic fields. In §3, we give results for the reionization process. In §4, we 
compute and discuss the global HI signal. In §5 and 6, the fluctuating component of the HI signal 
is derived and discussed. The detectability of the resultant signal is also discussed. Section 7 is 
reserved for summary and conclusions. In this paper we use the parameters of the FRW model as 
favoured by the recent WMAP results spatially flat FRW model with Q m = 0.3 and ft a = 0.7 
with n b h 2 = 0.022 and h = 0.7 & 



2. Magnetic fields in the post-recombination era 

Sethi and Subramanian (2005) [j|] analysed two possible effects of the presence of primordial 
magnetic fields in the post-recombination era: (a) magnetic fields can induce formation of early 
structures, (b) the dissipation of magnetic field energy owing to ambipolar diffusion and decaying 
turbulence can significantly alter the thermal and ionization history. 

The main difference between the standard ACDM models and the magnetic field-aided struc- 
ture formation is that, in the latter case, the magnetic field induces additional matter density fluctua- 
tions whose power spectrum is given by: P(k) oc k n for k < kj. Here A; j ~ 15(10~ 9 G/-Bo)Mpc _1 
is the magnetic field Jeans' wave number. The matter power spectrum spectral index n tends to 
one as the magnetic field spectral index approaches —3 (for details and discussion see e.g. [13, 8|j). 



Throughout this paper we use magnetic field power spectrum index —2.9, which gives n ~ 1.1. 

This results in more power at smaller scales as compared to the standard ACDM case. In 
Paper I [g], we showed that the la collapse of structures is possible at large redshifts; the collapse 
redshift of dark matter haloes is sensitive to the spectral index of the magnetic field power spectrum 
even though it is fairly independent of the magnetic field strength. It was also shown that the mass 
dispersion cr(R) oc 1/R 2 for scales above the Jeans' mass scales, which is a much sharper fall 
than the usual case. This means that most of the power is concentrated around scales close to 
magnetic Jeans' scales. These considerations also give strong indication that the only allowed 
tangled primordial field models are the ones for which n ~ — 3([j|, 13]). 



The magnetic field dissipation in the post-recombination era can also significantly alter the 
thermal Jeans' length scale, by heating the matter up to T m ~ 10 4 K, depending on the value of 
the magnetic field strength. Therefore, in this scenario, there is a play-off between early structure 
formation owing to enhanced power at small scales and suppression of power owing to magnetic 
field Jeans' scales and the enhanced thermal Jeans' scale (for details see [0]). Taking into account 
all of these effects, the interesting range of magnetic field strength is: 3 x 10~ 10 G < Bq < 
3 x 10~ 9 G. We study here the effects of primordial magnetic fields which are in the above range. 
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3. Reionization 



We model the reionization of the universe as described in Sethi (2005) [|i3|](for more details see 
[16]). In this model, an HII sphere is carved around collapsed haloes. Reionization proceeds as 
more sources are born and as the radius of the HII region around each source increases. It is 
completed when the HII regions coalesce. In this scenario, the ionized fraction at any redshift is 
given by: 

fion(z) = % [ dz ' J dM ^? (M ' Z ' )jR3(M ' Z ' Z/) (3 " 1} 

This expression remains valid until f lon ~ 1. Here the mass function dn/dM is computed using the 
Press-Schechter formalism. The radius of the HII region R is computed by following its evolution 
around a source with a given photon luminosity N 1 (in sec -1 ) which is assumed to grow linearly 
with the halo mass; we fix the fiducial value of the photon luminosity to be iV 7 (0) at the mass scale 
M = 5 x 10 7 M (see e.g. [15]). The radius of the evolving HII region is given by (e.g. Jl7| , [T|]): 

dR HR = (Nj ~ (47r/3)R s a B Cn 2 b x m ) 
dt (iV 7 + 4irR 2 x m n b ) 

Here a# is the case B recombination coefficient, n& is the background density of baryons (exclud- 
ing the helium atoms), and C is the clumping factor. We further assume the photon luminosity to 
be exponentially suppressed with time constant of one-tenth of the Hubble time at the redshift at 
which the source is born. The other parameter used in the evolution of the HII region around each 
source is the average clumping factor C; we assume C = 2 in the entire redshift range. 

As noted above, the mass function of objects in the magnetic field-aided collapse of structures 
is expected to have much larger number of objects close to the magnetic Jeans' scale as compared 
to the usual case, with this number falling sharply at larger scales. Tashiro & Sugiyama (2006a) 
[§] computed the reionization history in this case and showed that reionization is possible with 
star-formation efficiency markedly lower than in the usual case. We re-look at this issue within the 
framework of the semi-analytic approach described above, which explicitly takes into account the 
growth of HII bubbles around sources. 

We show our results in Figure 1. In the Figure, we plot reionization history for two values 
of magnetic field strength and the usual case with zero magnetic fields. We adopt iV 7 (0) = {4 x 
10 48 , 10 49 }sec -1 for the two magnetic field cases shown in the figure, respectively with Bq = 
3 x 10~ 9 G and Bq = 10~ 9 G. For the standard case with zero magnetic field we take iV 7 (0) = 
2.5 x 10 50 . All the reionization histories are normalized to the WMAP result r re i on = 0.1. 3 

Photon luminosity can be cast in terms of the star formation efficiency, f e g and escape fraction 
of hydrogen-ionizing photons / csc from star-forming haloes (for more details see [|l8|], p57). A 
photon luminosity of iV 7 (0) = 2 x 10 50 corresponds roughly to / eS c/eff = 0-01 if trie photons 
are emitted over a period one tenth of the Hubble time at z ~ 10, using Scalo IMF. This means 
that reionization in the standard zero field case requires both the star formation efficiency and 
escape fraction of roughly 10%. The product of these two parameters could be nearly two orders of 



3 We note that r rc i on is the integrated value between 30 > z > 0, to be relevant for the WMAP results. Owing 
to magnetic field dissipation, partial reionization can lead to a much larger 'reionization' optical depth at much higher 
redshifts. However, WMAP results are not sensitive to such early partial ionization (Paper I) 
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magnitude smaller if i?o — 3 x 10~ 9 G. Thus one needs a much lower efficiency of star formation 
for models with primordial magnetic fields. As both these parameters are highly uncertain our 
analysis doesn't allow us to constrain the magnetic field models from WMAP results alone (e.g. 
JT8|]). However, we note here that future meaningful constraints on the star formation efficiency 
and escape fraction could be used to rule/bear out magnetic field-induced reionization. 
We discuss briefly each of the case of magnetic field-induced reionization: 
flo ~ 5 x 10~ 10 G: In this case, the objects that cause reionization are in the mass range 
< a few 10 7 M Q or these are molecular-cooled haloes (for a review see e.g. [18]). There are many 
aspects to the nature of reionization caused by these objects that deserve further attention: (a) the 
temperature of the regions of the universe that are already (photo)ionized is ~ 10 4 K; this results 
in a thermal Jeans' mass too high to allow for baryonic collapse and star formation in molecular- 
cooled haloes in these regions. This means that as the reionization proceeds, in an increasingly 
larger fraction of the universe, baryonic collapse and star formation is suppressed for haloes with 
masses £ 10 8 M . Thus there is a quenching in the formation of new star forming haloes as the 
reionization proceeds. We take into account this process in computing the reionization history, (b) 
It has been argued that the UV light from the first objects can destroy the molecular hydrogen in 
the neighbouring haloes, thereby shutting off the process of further collapse and star formation in 
these haloes [[l 7 ^]. In Sethi, Nath, and Subramanian (2008) [14], we showed that in the presence of 
primordial magnetic fields, the changed ionized history results in a sharp increase in the formation 
of molecular hydrogen both in the IGM and in the collapsing halo. For magnetic field strength 
Bo > 5 x 10~ 10 G, the molecular hydrogen in the IGM might be sufficient to block the UV light 
from penetrating the IGM ([^, 14]). (c) Owing to the dissipation of the magnetic field energy 
in the post recombination era, the change in the thermal history of the universe can increase the 
thermal Jeans' mass. We take into account this change in our analysis. 

B > 10~ 9 G: In these cases, the collapsing haloes are in the mass range which allow both 
molecular and atomic cooling. This means that the impact of quenching, destruction of molecular 
hydrogen, and the increase in thermal Jeans' length discussed above is less pronounced [14]. 



4. The global HI signal 



The global HI signal is observable at the redshifted HI hyperfine transition. The signal can be 
seen against the CMBR, the only radio source at high redshifts. The observable deviation from the 
CMBR temperature at any frequency vq 



AT, 



1420/(1 + z) MHz is (e.g. [|2T 

Thi 



15]): 



CMBR 



(1 + Z[ 



CMBR 



(4.1) 



The spin temperature T s is given by [ 23 ] : 

T _ gcMBR + UcTk + y a T a 

i + y c + y a 



(4.2) 



Here Tr is the matter temperature and T a is the colour temperature of the radiation close to the 
Lyman-a line. In the context of the problem at hand, T a = Tk [24, 25]. y a determines the 
contribution of Lyman-a radiation in the determination of the spin temperature; y a ~ before the 
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onset of the epoch of reionization. r m is given by (e.g. [22, 15] and reference therein) 



... 'OHitf \ / 0.15 V/ 2 /T CM br\ A + *\ 1/a 



Here SIhi = pm/pc is the density parameter corresponding to the neutral fraction of hydrogen. 

The main difference between the standard zero field case and the magnetic field model is in the 
thermal history of the universe. The dissipation of tangled magnetic fields can significantly alter the 
thermal history of the universe (Paper I). The other difference is that the build-up of the Lyman-a 
radiation during the reionization era depends on the mass function of collapsed haloes, which in 
our case is computed using the magnetic field-induced matter power spectrum. 'Lyman-a' here 
refers to the continuum radiation between Lyman-a and Lyman-/? frequencies. The photons in this 
frequency range escape into the medium without getting absorbed by either the emitting source or 
the HI in the IGM until these photons redshift to a frequency very close to the Lyman-a resonant 
line frequency (for details e.g. fli5| ] and references therein). We model this effect, as in Sethi (2005) 



[15], by following the evolution mean specific intensity of the Lyman-a radiation: 



'"•W^ 5 ^/ mJ^^^ ( M,^ l + ^ (4.4, 

Here N ia is the Lyman-a photon luminosity of an object and (3 is the spectral index of the photons 
in the frequency range between Lyman-a and the Lyman-/5. Here 1 + z max = (1 + z)vp/v a and 
is determined from the knowledge that photons with an energy above the Lyman- j3 are absorbed 



locally (see e.g. [ ]26| , [27[]). Given the uncertainty about the spectrum of the ionizing sources, we 
assume = and N ia = ANj i.e. the 'Lyman-a' photon luminosity is a constant multiple of 
the luminosity of the Hydrogen-ionizing photons (for more details see the discussion following 
Eq. (10) in Sethi 2005 

We show in Figure 2 the global HI signal in the post-recombination era for A = 20, the ratio of 
Lyman-a to hydrogen-ionizing luminosity. The HI signal for one model with zero magnetic field 
is also shown; the zero magnetic field case shown corresponds to A = 20 with spectral index at 



energies exceeding Ly man-limit a = 2 (for further details see [15]), and is normalized to the three 
year WMAP results. The figure shows that unlike the usual case with zero magnetic fields, the HI 
signal in the magnetised universe is only observable in emission throughout the post-recombination 
era. This is mainly because the matter temperature doesn't falls below the CMBR temperature for 



the magnetic field models we consider (see also Ql 1|[) 



Different features of Figure 2 can be understood in terms of the altered thermal history in the 
presence of the magnetic field. The dissipation of the magnetic in the post-recombination era (for 
details see [0]) raises the matter temperature above the CMBR temperature. As the spin temperature 



(Eq. |4.2| ) is determined by matter temperature at high redshifts, the HI signal is observable in 
emission at high redshifts. This should be contrasted with the standard case in which the HI is 
observable only in absorption during the pre-reionization era (Figure 2). 

As the redshifts decreases, the dilution of the gas density drives the spin temperature towards 
CMBR temperature. The HI signal therefore approaches zero before the onset of the epoch of 
reionization. As the reionization epoch sets in, the spin temperature gets coupled to the Lyman-a 
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radiation, which results in the spin temperature being driven towards matter temperature again. In 
the standard case, depending on the details of the onset of reionization, the HI could be observable 
in absorption just as the reionization begins, as is seen in Figure 2 (see e.g. [15, 28]). However, in 
the magnetic field case, the signal is only observable in emission, as the matter temperature remains 
above the CMBR temperature. This part of the signal is potentially more promising, as it occurs at 
frequencies more accessible to the on-going and future experiments. The lack of detection of the 
HI in absorption therefore could signal the presence of magnetic fields. 

In both the standard as well the magnetic field case, the HI signal is observable in emission 
during most of the phase of reionization, with the signal strength almost independent of spin tem- 
perature as the limit T s » Tcmbr is reached. 



5. Fluctuating component of HI signal 



Here we analyse this signal in terms of two-point correlation function of the HI signal in real space. 



The correlation function can be expressed as (see e.g. [|29|]: 



with 



C(r 12 ,9) = (AT(n 1 ,r 1 )AT(h 2 ,r 2 )) = 

= To(ri)T (r2)[{^(n 1 ,r 1 )ip(n2,r 2 )) - (V>) 2 ] 

^ = x H (r)[l + <5(r)]. 



(5.1) 



(5.2) 



Here r\ 2 = |n — r2| and 6 is the angle between the line of sight h and the vector ri — r 2 
separating the two points, and brackets denote average over ri (the correlation function could be 
written equivalently as a function of the separation angle between the two points and the frequency 
difference u 2 — v\). 



[29, 30]) 



The expression in the square brackets in Eq. (5.1) can be further simplified to (for details see 

{iH>) - (V) 2 = txxfa(rvi, 0, z) + Cxx ~ x 2 H (5.3) 

Here £$8 = (£( r i)^( r 2)) is the two-point correlation function of the total density contrast, and 
ixx = (^h^i^h^)) is the correlation function owing to the inhomogeneities of the neutral 



fraction. For further detail on the assumptions made in deriving Eqs. fl5.1)-fl5.3) and the meaning 



of relevant variables see e.g. [30]. 



5.1 Density and ionization correlation functions 



The density correlation function, in the linear regime, can be expressed as (e.g. [ ]32| , [31| , [300) : 

is5ir X2 ,e, z) ~ £ ss (r 12 , 0, z) M + lp + X -(i 2 J + 



-K,vv(n2.2.:)(^ + ^ 2 ) P,(H) 



(5.4) 
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Here (3 ~ /b, b is the bias of HI distribution with respect to the matter distribution dominated 
by the cold dark matter. P2{0) is the Legendre function with 1 = 2, and 

& 5 (r,0,z) = ^fl J dkk 2 P HI (k)j (kr) (5.5) 

tis6(r,2,z) = -^tt J dkk 2 P HI (k)j 2 (kr). (5.6) 

Here D + (z) is the growing mode of density perturbations and PHi{k) is the HI power spectrum; we 
compute it as: Pni(k) = b 2 P(k), where P(k) is the underlying matter power spectrum. Through- 
out we use 6 = 1. 

As discussed above, we envisage the growth of reionization as expansion and coalescence of 
spherical ionized bubbles around each source. In this scenario, the two-point correlation function 
of the HI distribution is non-zero even in the absence of density perturbations. The scale of ionized 
inhomogeneities is determined by the scale of the ionized bubbles. 

If we assume all these spherical bubbles to be randomly distributed and non-overlapping (we 
will discuss this assumption below), then for bubble of size R, one can derive [HH: 



(5.7) 



(1 - PsameOR)) ' 

with 

2W = (l-x H )f(r,R). (5.8) 

and 

3r ( ^2 1 2 



f ( r ' R) = 1 -m{ R i2 r )> (5 - 9) 

for r < 2R and zero otherwise. Here p sa ,me is the probability that both points in the computation 
of the two point correlation function lie in the same ionized bubble. This formalism can be readily 
extended to bubbles of difference sizes if the bubbles are non-overlapping by adding contributions 
from bubbles of different sizes (see e.g. [[30|]): 

R 

Here g(R) is the distribution function of the radii of the HII regions, normalized such that its 
integral over all bubble sizes is the average ionized fraction: f g(R) = /i on - 

We show in Figure 3 the distribution function of bubble sizes g(R) at epochs when the average 
ionized fraction is fi on = 0.5 (left panel) and fi on = 0.15 (right panel) (at redshifts z ~ {10, 12.5}, 
respectively; see Figure 1), for models with magnetic field values of Bq = 5 x 10~ 10 G (solid 
curves) and Bq = 3 x 10~ 9 G (dashed curves). The dashed-dotted curves correspond to the zero 
magnetic field case. 

From Figure 3, we notice the following main differences between the magnetised and nonmag- 
netised models: (a) the sizes of bubbles are much smaller in the case of the magnetic field-induced 
reionization and (b) the distribution of bubbles radii is far broader in the standard zero magnetic 
field case. Both these differences are due to the difference in the matter power spectrum in the 
two cases. As already noted above, most of the haloes are formed close to the magnetic Jeans' 
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length in the magnetic field case. In the usual case, the power is distributed over a wider range of 
scales. Also, as pointed out above, the scale of HII regions is smaller in the magnetic field case 
because most haloes are close to l-a events and therefore the number of collapsed haloes is larger 
and the iV 7 is smaller, as compared to the usual case. An important result of this paper is that the 
size distribution of HII regions is almost independent of the magnetic field strength as also seen in 
Figure 3. Even though the mass of the collapsed haloes are larger with increasing magnetic field 
strength, their luminosities have to be appropriately lowered to normalize to WMAP results, as 
already noted in a previous section. The results shown in Figure 3 for the magnetic field models 
are therefore generic and we note here that we did not find any exception to these features in a wide 
variety of models we studied. 

5.2 Effect of Large-scale density field on ionization inhomogeneities 

The results of Figure 3 implicitly assume that the positions of ionizing sources are uncorrected. 
This however is generally a poor assumption. In the usual case, it has been shown that the scale 
of bubbles is generally much larger than the assumptions of Figure 3, and the scale of ionization 
bubble also evolves very rapidly (e.g. [|33|, |34|]). This is also known to be in agreement with 
simulations (e.g. [37, |3f]]). In the usual AC DM case, the local density power spectral index 



in the range of scales of interest to the reionization generally lies in the range —3 to —1, or the 
density field is strongly correlated for a large range of scales. This means that the impact of the 
large scale density field on the HII region size is pronounced. Also, the first objects to form are rare 
objects (~ 2.5a), which means the centers of ionizing sources are more strongly correlated than 
the underlying density field. 

In the magnetic field case, the density power spectral index is ~ 1 above the magnetic Jeans' 
mass scale, or the density perturbations are strongly suppressed on scales larger than this scale. 
Also the first objects to collapse are close to la events which means the centers of the ionizing 
sources are less correlated. This means that the effect of the large scale density field in this case is 
likely to be less pronounced. However, we still need to take it into account to assess its impact. We 



follow here the prescription of Furlanetto, Zaldarriaga, and Hernquist (2004) Q34{]. 



Furlanetto et al. (2004) [ 34 ] showed that the comoving number density of the HII regions with 



masses in the range M and M + dM can be expressed as: 

dn 2 Oh dlna Bn /•„•>, * , , o „ \ 

K«7 exp W M '«' (5 " n) 



M dM 



Here B(M,z) = B + B x a 2 with B = 5 c (z) - V2K(()a m - m and B x = K{Q/{^2a~); 
K{C) = erf _1 (l — and cr m i n corresponds to the mass dispersion at the size of the smallest 
halo that can cause reionization. We assume this mass to be M ~ 5 x 10 7 M throughout. The 
parameter Q is the ratio between the ionized fraction and the collapsed fraction at any given redshift. 
As C tends to one, the mass function of the HII regions approaches the standard Press-Schechter 
mass function. This is the limit in which the HII regions enclose a single ionizing sources, and 
therefore is expected to approximate the results shown in Figure 3. However, ( is generally much 
greater than one, which means the HII regions enclose multiple ionizing sources, and their sizes 
evolve far more rapidly than the prescription of Figure 3. 
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The size distribution of HII bubbles is given by the single parameter ( in this formulation. 
This parameter would generally depend on various physical parameters e.g. photon luminosity 
of collapsed haloes, clumping factor, etc, and in general will evolve depending on the history of 
reionization. In other words, the computation of HII regions sizes would require a physical model 
of the evolution of the ionized fraction. In this paper, we compute this parameter as the ratio of 
ionized fraction to the collapsed fraction at any redshift from the the semi-analytic model described 
in section 2, for the ionization history shown in Figure 1. 

In Figure 4, we show the distribution of HII region sizes for two values of fractional ionization, 
/i on and three values of magnetic field strength. The values of £, in increasing order of the strength 
of the magnetic field, are: ~ {10, 14}, {9, 12}, {4, 7}, with each pair of values corresponding to 
the ionized fraction /j on ~ {0.15, 0.5}. These results should be compared to Figure 3. We notice 
that HII regions sizes are larger and evolve more strongly. However, as compared to the usual case 
(e.g. Figure 4 of [34]), the evolution is not as strong for the reasons discussed earlier in the section. 

We use the results of Figure 4 to compute the fluctuating component of the HI signal owing to 
ionization inhomogeneities. 



6. Results: Fluctuating component 

We show the two-point correlation function of the HI fluctuation in Figure 5 for several magnetic 



field models. All the curves are for 9 = ir/2 (Eq. (5A)), that is they basically focus on the angular 
correlation function. We also break the contribution owing to the density and ionization inhomo- 
geneities (the former corresponds to the first term in Eq. ( |5J| ) and the latter to the next two terms 
in the equation). 

The first three panels of Figure 5 (counterclockwise from bottom) correspond to signals which 
arise for magnetic field values: {Bq = 5 x 10~ 10 , 10~ 9 , 3 x 10 -9 } G, respectively. The dashed 
curves give the contribution of the first term in Eq. (p3|), that is £ xx £f s , where is the correlation 
function of the purely magnetically induced density perturbations. The dot-dashed curves show the 
contribution of ionization inhomogeneities to the HI correlation function (the £ xx — x| term). The 
solid curves in these three panels give the sum of both these contributions. 

Note that even in the models with primordial fields, there is another contribution to HI fluctua- 
tions (and correlations) due to the inflationary induced density perturbations, say referred to as 
It should further be pointed out that if the reionization is caused by primordial magnetic field, the 
photon luminosity of haloes is too small (see discussion on Figure 1 above) for the haloes formed 
owing to inflationary density perturbation to cause significant contribution. The number of these 
haloes is too small to be important in the process of reionization. Therefore, the only impact of 
the inflationary perturbation in the presence of primordial magnetic fields is to add a contribution 
owing to density perturbations. In the fourth panel (top left panel of Figure 5), we show the contri- 
bution ( ^xx^sf) as a dotted line. The solid line in the fourth panel shows the result of adding this 
contribution as well to the magnetically induced signals with Bq = 3 x 10~ 9 G. Thus it shows the 
predicted total HI correlation function from all sources of density and ionization fluctuations for 
the case Bq = 3 x 10~ 9 G. 

The results shown in Figure 5 can be summarized as follows: 
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1. The typical scale of ionization inhomogeneities is ^ 1 Mpc for the magnetic field models. 
This results follows from the nature of distribution function of HII regions, as shown in 
Figure 4. 

2. Ionization inhomogeneities make a comparable contribution to HI fluctuation signals as that 
arising from density perturbations. For smaller values of magnetic fields, the ionization 
inhomogeneities can also provide the dominant contribution of the signal. 

3. The HI correlation function arising from density fluctuations show characteristic oscillations, 
with a scale length which increases with increasing value of the magnetic field strength. This 
is expected as the dominant contribution to the signal comes from scales close to the mag- 
netic Jeans' length. Detecting such oscillations could give an indication of the influence of 
primordial magnetic fields and also help in determining the field strength. Also the change in 
sign of the HI correlation at scales above the scale of ionization inhomogeneities (r ^ 1 Mpc 
could provide additional information. 

4. In Panel 4 of Figure 5 (counter clockwise), we show the total observable signal in the pres- 
ence of a magnetic field of strength 3 x 10~ 9 G. As noted above if the magnetic field-induced 
structure formation is responsible for reionization, £ xx is totally determined by the primor- 
dial field strength, and the only contribution of inflationary generated density perturbations 
is to add a term of the form ^xx^sf to tne HI correlations. 

The main difference between the usual AC DM case and the magnetic field scenario, as also 
alluded above, is in the scale of the signal. Both the density coherence scale (as seen in the fourth 
panel of Figure 5) and the scale of ionization inhomogeneities are larger in the ACDM case. 

The process of reionization can be considerably more complex as compared to the results 
shown in Figure 1 . The ionizing sources can have varying spectral indices, star-formation efficien- 
cies and star formation histories, which could have a non-trivial impact on the photon luminosity 
in the process of reionization. We have studied a suite of ionization histories in addition to the 
one shown in Figure 1 and find that the results in Figure 3 and 4 are fairly robust. This behaviour 
arises due to the fact that the magnetic field scenario is completely dominated by scales around the 
magnetic Jeans' length. This scale determines both the density inhomogeneities as well as the sizes 
of HII regions that form. 

6.1 Detectability of the signal 

As noted above, the main difference between the standard zero field and the magnetic field induced 
reionization is generally the much shorter correlation length scale of the fluctuating HI signal in the 
latter case. Many radio interferometers currently being built to observe the epoch of reionization 
have typical angular resolution of a few arc-minutes (1' ~ 2 Mpc for the cosmological parameters 
of interest), which corresponds to a linear scale ~ 5-10 Mpc. 4 . These experiments are therefore 
not sensitive to the scales at which magnetic field effects are likely to dominate, except for the 
largest values of Bq. MWA is expected to have a primary beam of 20° and synthesized beam of 

4 These experiments have greater resolution in the frequency domain, e.g. MWA can have 4000 channels over 32 
MHz band-width, which corresponds to length resolution of ~ 15 kpc 
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4' with 500 'tiles'. Assuming full UV coverage of the primary beam, the expected error on the 
estimated two-point correlation function is ~ 3 x 10~ 7 K 2 , for an integration time of 10 6 sec (see 
e.g. [30]). This might allow a detection of the signal, as shown in Figure 5, with magnetic field 
of strength Bq ~ 3 x 10~ 9 . However, we note here that an indirect indication of the presence 
of magnetic field might be possible. If the reionization was indeed caused by a magnetic field of 
strength Bq ~ 10~ 9 G, the reionization process would appear to be almost 'homogeneous' at the 
angular resolution probed by these instruments. 

Future radio interferometer SKA has the angular resolution and the frequency coverage to 
directly detect the redshifted HI fluctuation owing to magnetic fields effects, for the entire range of 
magnetic field strengths we consider here. SKA is projected to have a point-source sensitivity of 
~ 400 /j, Jy in one-minute integration (in the continuum mode) in the frequency range of interest 
(70 MHz and 300 MHz). Assuming a synthesized beam in the range 15"-30", relevant for our 
study, and a simultaneous sky coverage of 25° (the simultaneous sky coverage of SKA is planned 
to be much larger, which will increase the sensitivity further), the sky brightness sensitivity for 
the statistical detection (in the line mode) is in the range 10~ 6 K 2 -2 x 10~ 7 K 2 for one month of 
integration. This suggests that a 5-a detection of the signal shown in Figure 4 should be possible 
with less than one week of integration with SKA. 



7. Summary and Conclusions 



Tangled primordial magnetic field can have interesting consequences for cosmology. The presence 
of these fields leave detectable signatures in the CMBR temperature and polarization anisotropics 
(pi, M, pf|, 0],©, [@,[0]J0],@],@,ra,@il). The statistics of this signature will 



also have distinctive non-Gaussian features [51, 52, 53]. From these considerations, one can obtain 



an upper limit on the magnetic fields strength Bq ^ 4 x 10 9 G [43], for models with magnetic 



field power spectrum index n ~ — 3, which are the only models that are consistent with current 
state of observations [54, 80. 

Sethi & Subramanian (Paper I) considered many aspects of these fields in the post-recombination 
universe. In particular, they showed that magnetic field-induced structure formation might lead to 
early reionization and the dissipation of these fields can substantially alter the thermal and ioniza- 
tion history of the universe. 

Magnetic field-induced reionization and its possible impact on the detectable HI signal has 
been considered by several authors [[[(J, O]. In this paper, we study these aspects further within the 
framework of semi-analytic models which allow us to compute and add in the impact of reionization 
inhomogeneities on the HI signal. 

Our main findings can be summarized as: 

1 . Owing to magnetic field dissipation in the post-recombination, the HI signal is not expected 
to show any absorption features unlike the standard case (Figure 2) (see also [|P1|]). 

2. The matter power spectrum owing to magnetic-field induced structure formation is peaked 
around the magnetic Jeans' scale. This scale is imprinted on the HI signal as typical oscilla- 
tions in the HI angular correlation function (Figure 5) 



- 11 - 



3. Figure 5 also captures the impact of ionization inhomogeneities. The scale of the ionization 
inhomogeneities is < 1 Mpc, for a wide range of magnetic field strengths (Figure 4). 

4. We show that the presently-operational and upcoming radio interferometers are sensitive to 
the magnetic field strength 5o~3xlCP 9 G. However, the future SKA will be able to detect 
the HI signal for the entire range of magnetic field strength: 5 x 1CT 10 ^ Bq ^ 3 x 1CT 9 G. 
A detection of a characteristic scale and nature of oscillations of the HI correlation function 
could help measure the magnetic field at these scales directly. 

The upcoming SKA will also have the capability of direct detection of primordial tangled 
magnetic field by measuring the Faraday rotation of ~ 10 7 radio sources. 

We note here that the main advantage of a potential detection of the HI signal is that it is 
sensitive to far smaller values of magnetic field strength as compared to the other probes. 
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Figure 1: The ionization history of the universe is shown for two magnetic field models along with the usual 
case without magnetic field. The dashed and the dot-dashed curve correspond to the magnetic field strength 
£?o = {10~ 9 , 3 x 10~ 9 } G, respectively. The solid curve correspond to the case without magnetic field. For 
details about other parameters see the text. 
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Figure 2: The Global HI signal is shown for two values of magnetic field strengths. The solid and dashed 
curves correspond to the magnetic field strength B = {5 x 10~ 10 , 10~ 9 } G, respectively. The dot-dashed 
curve corresponds to HI signal for one possible scenario in the zero magnetic field case (see text for details). 
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Figure 3: The ionized fraction is shown as function of radius of the HII region. The area under each curves 
is the ionized fraction. The dot-dashed curve correspond to the case without magnetic field. The solid and 
dashed curves correspond to the magnetic field strength B Q = {5 x 1CT 10 , 3 x 10~ 9 } G, respectively. In the 
left panel the ionized fraction is /; on ~ 0.5 and /i on ~ 0.15 for the right panel, which correspond roughly 
to the redshifts 10 and 12.5, respectively, as discussed in the text. 
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Figure 4: The ionized fraction is shown as function of radius of the HII region. The solid, dashed, and 
dot-dashed curves correspond to {B = 5 x 1CP 10 , 1CP 9 ,3 x 10~ 9 }G, respectively. The lower and the 
upper curve for each line style correspond to /i on — 0.15 and /i on — 0.5, respectively. ( values for these 
curves are: ~ {10, 14}, {9, 12}, {4, 7}, with each pair of values corresponding to the two cases of ionized 
fraction. 
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Figure 5: The first three panels (counterclockwise from bottom) correspond to the signal for magnetic field 
values: {B = 5 x 10~ 10 ,10~ 9 ,3 x 10~ 9 } G, respectively; The signal corresponds to / ion ~ 0.5 (Left 
panel of Figure 3) at z ~ 10. The solid, dashed, and dot-dashed curves correspond to the absolute value of 
the total signal, the signal owing to density fluctuations, and the signal from ionization inhomogeneities. In 
the forth panel, we show the total observable signal (solid line) for B = 3 x 10~ 9 G, which includes the 
contribution from density fluctuations in the standard ACDM case (dotted line). 
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